Incommensuration in quantum antiferromagnetic chain 
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A dimerized quantum Heisenberg or XY antiferromagnetic chain has a gap in the spectrum. 
We show that a weak incommensurate modulation around a dimerized chain produces a zero tem- 
perature quantum critical point. As the incommensuration wavelength is varied, there is a transi- 
tion to a modulated gapless state. The critical behaviour is in the universality class of the classi- 
cal commensurate-incommensurate (Pokrovsky-Talapov) transition. An analogous metal-insulator 
transition can also take place for an incommensurate chain. 
The one dimensional 5 = 1/2 quantum antiferromagnetic Heisenberg model is a prototypical example of quantum 
i-G ■ fluctuations destroying the classical ground state Such a uniform chain, with gapless excitations, can undergo a 

| spin-Peierls (SP) transition to a lattice dimerized state via the interplay of spins and lattice distortion, going over to 
an incommensurate phase in large magnetic fields Jz|— p~0[ . The SP state has a gap in the spectrum. However, specific 
heat measurements in the incommensurate phase of the inorganic insulator CuGeC>3 Jl(| hint Jll],[l^] towards a gapless 
phase. This, then, raises a few basic questions: Can a weak incommensuration around a dimerized Heisenberg (or 
XY) chain support gapless excitations, eventhough the dimerized case does not? If it is true, then is there a quantum 



^ . critical point (QCP) separating a gapless state and a gappy state (a state with a gap)? What would be the nature of 



the quantum critical point? 

In order to address these issues for the incommensurate phase, and to gain an understanding of the underlying phe- 
nomena, one may consider a simpler situation of a static case where the coupling constant follows the incommensurate 



| modulation of the lattice. We therefore consider a model hamiltonian 

C ' N 

o 
o 



H = H ^(S}Sr +1 + S 7Sf +1 ) +hJ2s*, (1) 
j'=i L J i 

where Sf = SJ ± iS!j are the spin raising and lowering operators at site j, and h is the magnetic field in the z- 
direction. The case of interest in this paper is Jj = J + <5cos(7r + q)j, with 6 <C J > 0. This corresponds to a weak 
incommensuration around the dimer phase (q = 0). For S = 0, the system is gapless (uniform chain). For q = 0, a 
| gap develops in zero field for 5^0 but a gapless phase is recovered above a critical magnetic field h c jl3| . Our aim 
OO ■ is to get the phases in the (q, h) plane for 5 <C 1, based on low energy excitations. 

We are considering the XY model, because, so far as the gap is concerned the difference between XY and Heisenberg 
chain is not crucial (see Eq. (|j)). Admittedly, for the SP problem, the spin lattice interaction and the magnetic field 
play important roles Jlif and correlate Jlo| q and h. A real system would follow a particular trajectory q = q(h) in 
the (q,h) phase diagram. We therefore take q and h as independent parameters. 

Eq. (Q) also describes, via a Jordon Wigner transformation spinless fermions on a one dimensional lattice 
- a description we use in this paper. By defining Cj — K{j)Sj as the fermionic annihilation operator, where 

K(j) = exp(«7r^^~ 1 1 S+S~) is the nonlocal kink operator, the spin hamiltonian of Eq. [j]can be written as 

o 

o ■ 



H = \ E J i(C]C j+1 + Cj +1 C 3 ) -hJ^Cfo. (2) 



' Note that the fermions have incommensurate hopping rates, and h acts as a chemical potential. 

A QCP is a phase transition point with a diverging length scale induced by a change of a parameter of the 
hamiltonian at zero temperature. Since quantum fluctuation is responsible for such criticality, the dynamic exponent 
z, determining the scaling of time and space, becomes more important than in thermal critical points (r ~ £ z with 
t and £ as characteristic diverging time and length scales). In the spin chain context, a wellknown QCP is the 
dimerization point S = for q = h = for Eq. (|l|) or (J2J) . This critical point (separating a gapless state and a gappy 
state) corresponds to free fermions with z — 1 and a correlation length diverging as £ = The spin-spin correlation 
function C(r,5) = (S*S? +r ) has a scaling behaviour C(r,S) = r^^f^rS) with n = 2 Jl7| ). For the QCP of concern in 
this paper, we start from the dimerized gappy state and change the incommesuration q keeping 8 fixed. 

A QCP for q — at S = ensures the existence of a continuum limit of the lattice problem. We use this continuum 
limit to study the effect of q on the phase diagram, and in a renormalzation group framework only the relevant terms 



f 



need to be kept. To go to the continuum limit, we adopt the technique of bosonization for the low energy excitations 
around the fermi points |]l8| , |l6| , ^9| . Also since wave number is nolonger a good quantum number, we work in the real 
space. The low energy excitations of Eq. (0) across the fermi surface on the two branches can be described by L (left) 
or R (right) moving particles with linear spectrum. This gives z = 1. In the continuum limit, the Hamiltonian is 
written in terms of the (bosonic) phase variables. The free bosonic theory corresponding to Jj = constant, h = turns 
out to be just the harmonic lattice hamiltonian, which forms the basis for the subsequent analysis of the remaining 
terms such as dimerization, incommensuration around the dimerized case and the magnetic field. The relevance or 
irrelevance of various terms come from their scaling behaviour on a long length and time scale. In the renormalization 
group approach, the scaling behaviour is obtained by integrating out the short distance fluctuations in both directions 
and incorporating their effects in the parameters of the problem. This is implemented by studying the Euclidean 
version of the problem (i.e., imaginary time, t — > it). 

For low energy excitations, we define new operators Oj = (i)~ J Cj to eliminate the fast variation at the Fermi 
vector. From these, right and left movers are defined by the chiral transformation R(j) — (a 2 j — a,2j~i)/\ // 2 and 
L(j) = (a 2 j + &2j_i)/\/2 so that the continuum version of the Hamiltonian can be written as 



H = i 



J dxJ[L\x)d x L{x) - R\x)d x R{x)] - (3a) 



- / dx dJ(x){R j '(x)R(x) - L f (x)L(x)) + (3b) 
% -a I dxd 2 J(x)(R\x)L(x) - tf(x)R{x)) (3c) 



where for simplicity the modulated part of the coupling has been ignored in Eq. 3a, This is justified by the modified 
Harris criterion of Ref. |2l| since incommensuration 5cos(2kf + q)x is bounded. The notation dJ{x) and d 2 J(x) 
denote the appropriate continuum limit of the discrete difference of the modulation of the couplings. 
The bosonic operators are obtained from the phases ||^|| of R(x) and L(x), 

R{x) ^ —^^=exp(ik f x ~ ^(x)), (4) 

L(x) = —j^=exp(-ik f x- (f} 2 {x)) 7 (5) 

with kf = 7r/2, so that 9(x) = i(4>i(x) + ^(x)), and the conjugate momentum p are related to the density and current 
respectively as 89/ dx = R^R + L^L and p = R^R — L^L. With this choice, the free part of the bosonic hamiltonian 
(analogous to the first term, Eq. pal) is the standard harmonic chain hamiltonian 



H = ±- 
2vr 
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j dx 
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\dx) 



(6) 



in units of h = 1 and spin wave velocity = 1. Here, K = 4 for the XY model. For the full isotropic Heisenberg 



hamiltonian or interacting fermions in the fcrmionic language, K = 2. The second term, Eq. 3b, is analogous to 
a current type term which is absent if one goes over to the dimer limit (recovering translational symmetry) and is 
expected in the incommensurate case on symmetry grounds (or lack of it). In the bosonic variables, the current term in 



Eq. 3b is like 5 q sin Qx p(x), where Q — 2kf + q, and it vanishes in the limit q — > 0. However because of the oscillatory 
coefficient with average zero over a period, we ignore this term in the present analysis. The third term ( |3cl ) shows that 
a spatially varying J gives rise to a umklapp scattering process which apparently violates momentum conservation 
by 2k f. This is a special term needed mainly because of the underlying lattice in the problem. Considering only slow 
variations in the bosonic field 0, the umklapp term can be wr itten |H as » ^ / dxQ 2 cos(9(x) - [2k f - Q]x), by 
shifting 9 — * 9 — n/2. The bosonized version of the incommensurate Hamiltonian is therefore (prime denoting spatial 
derivative) 



cos(9(x) + qx) (7) 



The dimer limit is restored by taking q — and in this limit with kf — -k/2 the umklapp term agrees with Nakano 
and Fukuyama [^2| who obtained a similar term in the presence of a constant bond alteration. The above bosonic 
hamiltonian has resemblance with the hamiltonian that appears in theory of incommensurate crystals (231] . In those 



2 



problems q plays the same role of incommensuration vector. Such a hamiltonian also occurs in the Frenkel-Kontorowa 
model p4| of a harmonic chain in an external cosine potential. This similarity shows that the solitons expected in 
such cosine potentials are also important in the spin chain problem. 

To study the relevance or irrelevance of the sine Gordon term, we implement a renormalization group analysis, 
well documented in Ref. p5| , p3| ]. The basic steps involve obtaining a functional integral J e lS ^d9 where S(9) is 
a dimensionless action. By going over to imaginary time t = it the integral is written as a partition function of 
a classical two-dimensional problem. A further shift of the field variable 6(x,t) = 9(x,t) — qx, gives the classical 
hamiltonian [Eq] as 



H = / drdx 



le 2 + ±(e'-q) 2 -Acost 



(8) 



where A = Stt 2 K/A. 

The RG procedure |25| for q = involves decomposing the cosine term into fast and slowly varying components of 9 
for A -C 1 and then averaging out the fast varying component with respect to its Gaussian distribution. Absorbing the 
contribution from the averaged fast component in A we can have the renormalization of the A as Ar — A 1 /' 1-73 / 87r ), 
where (3 2 = nK. For /3 < V&n, the sine Gordon term is relevant and yields a massive theory. The Heisenberg or the XY 
model belongs to this category and, as expected, any dimerization produces a gap in the spectrum, with the gap scaling 
as An, and therefore scaling as S 8 ^ 8 ^ K \ With the incommensurate term, the system can disregard the potential if 
q9' is comparable to the cosine potential energy. This gives the phase transition point Q as q c ~ < 5 4 /( 8 ^ Ar ). For the 
XY model, the transition to the gapless phase is in the same universality class as the commensurate-incommensurate 
transition in two-dimensions p3Lp7|. For this universality class, the relativistic invariance leading to z = 1 is lost. 
The spatial length scale exponent is v = 1/2 with z — 2 p9| . p8| . 

The results imply that the gap at the fermi level persists for small q incommensurations around the dimerized 
phase. Isolated states might well appear in the gap but they do not destroy the gap. The spins basically follow a 
dimerized lattice. For larger incommensuration, i.e. for q > q c , the spins follow the dimerized lattice over a finite 
length scale separated by defects or solitons. Such (spin 1/2) solitons do exist for the dimerized chain [p2|. Though 
the solitons cost energy, a many soliton state for large q would be favorable compared to, say, following the dimerized 
lattice over the whole length (a no soliton state) because of the energy gain through the "q" term in Eq. [|. In this 
phase with a finite density of solitons, translational invariance is recovered because these solitons are not bound to 
the underlying lattice, and one gets back a massless mode. The density p of these solitons is given by (9'), and this 
density can act as the order parameter for the transition. In the gappy phase, p is zero (no soliton or domain wall), 
while p is nonzero in the gapless phase. The dependence of p on the deviation from the critical point is obtained by 
matching the chemical potential to the relativistic fermi energy. This leads to a density ~ (q — q c )@ with (3 — 1/2. 
In one dimension, a length scale I can also be defined from p as the average separation of the solitons. This length I 
diverges (I ~ (q — <7c)") with an exponent v = j3 = 1/2 as the critical point is reached on the gapless side. Because 
of the high energy involved in the soliton formation, there will be no critical divergences in the gappy state |2^,^,^] . 
The gapless phase is not identical to the free fermion phase, mainly because of the existence of a length scale I for the 
average domain size within which the system follows the dimerized lattice. E.g., the equal time spin spin correlation 
function C(r, q) would have an algebraic decay as for free fermions or the dimerization QCP with rj = 2, but for the 
incommensurate QCP, there will be an additional oscillatory factor cos(r//) |^,^]. Such oscillatory factors would be 
detectable in scattering experiments. 

For the spin system, the phase transition induced by the incommensuration q signals the formation of a band at the 
fermi surface (in zero field) through the wandering of the solitons which in one dimension also act like noninteracting 
fermions p8| , p7| . The low lying spectrum therefore becomes gapless. A simple dimensional analysis suggests that 
the width of the band formed around the fermi level is w ~ l/l 2 . We, therefore, expect the bandwidth to vanish 
as w ~ |qr — q c | 2 ^,i.e., linearly on the gapless side. According to the bosonization rules the local magnetization is 
proportional to (9'} so that the soliton density also determines the magnetization in the phase. Therefore the gappy 
phase will be nonmagnetic but the gapless phase is magnetic and the magnetization vanishies with the Pokrovsky- 
Talapov exponent (3 — 1/2. 

So far we have considered the zero field case. The magnetic field acts as the chemical potential of the fermions, 
as seen in Eq. ^, and so long as the fermi surface is in a band, the equivalent bosonized hamiltonian will be of the 
form Eq. |^, with a renormalized K and the spin wave velocity ( taken to be unity). Since h couples to 89 /dx, such 
a magnetic field for a dimerized chain would have an equivalent bosonized hamiltonian as Eq. || with h replacing q. 
The critical behaviour as the fermi surface reaches the boundary of the gap will therefore be similar to what we have 
already studied. This has explicitly been shown in Ref. If the band formed by the solitons lies entirely in the 

original gap, then by shifting the fermi level the system can go from a gapless to a gappy state and then again to a 
gapless phase. 
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Our analysis though aimed at the spin problem is equally valid for hopping spinless fermions on a one dimensional 
incommensurate lattice i.e. on a lattice with incommensurate hopping rates, as given by Eq. ^. We can therefore 
conclude that for a half filled lattice, there will be an insulator-metal transition as the incommensuration wavelength 
is varied. One dimensional incommensurate systems can therefore be classified as metals or insulators based on the 
incommensuration. It might be possible to observe such a metal insulator transition in an incommensurate crystal 
(incommensuration around dimerized lattice) by changing say temperature or other external parameters that control 
the incommensuration of the lattice. 

As emphasized at the beginning, our analysis is tied to the dimerized gappy phase, the existence of a continuum 
limit, and the perturbative renormalization group (i.e. S <C 1). From the similarity with the Frenkel-Kontorowa 
model in the continuum limit, Eq. (0). and the possibility of a gap in the original Hamiltonian of Eq. [l] whenever the 
wavevector Q = n + q is a rational fraction of the fermi wavevector, it is expected to have stable lock-in phases [p3| 24 
around certain rational fractions. Our procedure would then yield similar critical behaviour around each such lock-in 
gappy phase. To get the full phase diagram (and the possibility of a Devil's staircase [^3|,|4|) one needs to study the 
original lattice model as the strength of the incommensuration S is increased. Also strong incommensuration may 
completely destroy the band structure yielding point spectrum as known for quasiperiodic systems [fjfjll . These remain 
to be studied. 

To summarize, we have studied the effect of incommensuration in a quantum antiferromagnetic Heisenberg or 
XY chain and found a quantum critical behaviour induced by incommensurability q around dimerization. Using 
renormalization group results we find that, in zero magnetic field, with the increase of the lattice incommensurability 
the massive theory reduces to a massless theory through a continuous transition at zero temperature. This indicates 
a transition from a nonmagnetic gappy (as in the dimerized case) to a gapless magnetic phase. The critical behaviour 
of this transition is found to be in the same universality class as the two dimensional commensurate-incommensurate 
transition, with all singular features appearing only on the incommensurate gapless state. Correlation functions will 
have characteristic oscillatory factors that distinguish the incommensurate phase from the gapless free fermion phase. 
This prediction could be tested in synthetic one dimensional magnetic chains. We have also pointed out that our 
results are equally valid for spinless fermions on special types of incommensurate crystals. Therefore the prediction 
and critical nature of the metal-insulator transition induced by incommensuration around a dimerized lattice can be 
tested in properly fabricated incommensurate hetero junctions |31J| . 
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